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Abstract - We study the influence of disorder on the topological transition from a two- 
dimensional Dirac semi-metal to an insulating state. This transition is described as a continuous 
merging of two Dirac points leading to a semi-Dirac spectrum at the critical point. The latter is 
characterized by a dispersion relation linear in one direction and quadratic in the orthogonal one. 
Using the self-consistent Born approximation and renormalization group we calculate the density 
of states above, below and in the vicinity of the transition in the presence of different types of dis- 
order. Beyond the expected disorder smearing of the transition we find an intermediate disordered 
semi-Dirac phase. On one side this phase is separated from the insulating state by a continuous 
transition while on the other side it evolves through a crossover to the disordered Dirac phase. 



Introduction. — In the recent years there has been 
a growing interest for materials whose low energy degrees 
of freedom are described by two dimensional (2D) Dirac 
fermions. For instance, 2D Dirac points have been fomid 
in the energy spectrum of graphene [HH], 3D topological 
insulators [2111] such as the materials in the 612863 fam- 
ily [5] and strained HgTc [B], quasi-2D organic conductor 
a-(BEDT-TTF)2 I3 under pressure [THTU] and very tun- 
able molecular system assembled by atomic manipulation 
of carbon monoxide molecules on a copper surface [111 . 
Each Dirac point carries a topological charge which is ro- 
bust against perturbations [12] . This topological charge is 
related to the chirality of the electronic states, but can also 
be viewed as a Berry monopole. The Berry phase associ- 
ated with any closed path in reciprocal space is determined 
by the total topological charge of the Dirac points enclosed 
by the path. Moreover, in the case of two-dimensional 
models, excluding the surface states, according to the 
fcrmion doubling theorem the Dirac cones emerge by pairs 
of opposite topological charges. In this topological 
transition towards an insulator is possible which consists 
of a merging of the Dirac points of opposite topological 
charges [T3] • Moving the Dirac points in the Brillouin zone 
up to their fusion requires perturbations of the underlying 
microscopic Hamiltonian. This can be achieved, for exam- 



ple, by mechanical stretching or depositing adatoms that 
could give rise to a hopping anisotropy or enhancement 
of higher-order hopping amplitudes. Exactly at the topo- 
logical transition triggered by the merging of two Dirac 
points, a very peculiar semi-Dirac dispersion relation oc- 
curs, which is linear in one direction while being quadratic 
in the orthogonal one [T5HT5] . Recently such a merging of 
tunable Dirac points has been realized using a degenerate 
Fermi gas trapped in a 2D honeycomb optical lattice [T6] 
and in microwave experiments |17j . Moving and merging 
of two Dirac points in the organic conductor a-(BEDT- 
TTF)2 I3 has also been discussed in fTSllTO]. 

For 2D Dirac semi-metals, both the density of states 
(DOS) near the Dirac point and the transport properties 
may be strongly affected by the presence of disorder. Its 
effect crucially depends on properties of disorder such as 
preserved symmetries or the range of correlations. This 
has been intensively studied using different methods dur- 
ing last two decades [20lI26j . On the contrary, much less 
is known on the disorder effects close to and at the topo- 
logical semi-metal-insulator transition. Very recently, the 
Friedel oscillations induced in the semi-Dirac fermions by 
a single localized impurity have been studied in [57] . In the 
present paper we address the effect of multiple randomly 
distributed impurities of different types on the DOS in 
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Fig. 1: The energy spectrum of the clean system: The left 
panel corresponds to A < 0, the right panel to A > 0. The 
center panel is the semi-Dirac point A = 0. 



the vicinity of the topological transition. We focus on the 
DOS at low energy, near the original Dirac points, while 
the different problem of diffusion of semi-Dirac fermions 
at higher Fermi energy will be presented elsewhere [28] . 

Model. — We consider a 2D semi-metal with valence 
and conduction bands touching each other at two Dirac 
points. We assume that the Dirac points can evolve in 
reciprocal space upon variation of the parameter A which 
drives the topological merging transition: for the criti- 
cal value A = the two points merge as shown in fig- 
ure [1] The effective Hamiltonian describing the low en- 
ergy physics in the vicinity of the transition can be de- 
rived from the tight-binding model expanding around the 
merging point. It is expected to be general and reads [13] 



Ho 



A 



2m 



iCOydy. 



(1) 



At the critical point A = 0, Hamiltonian ([T]) describes the 
semi-Dirac fermions with a quadratic dispersion relation 
parametrized by the mass m in the x-direction, and linear 
dispersion with velocity c in the y-dircction. We assume 
that m > and set ft = 1, while ctq = I and the three cr^, 
/X = x,y,z are the Pauli matrices. The diagonalization 
of Hamiltonian ([!]) via Fourier transform gives the energy 
spectrum e = ±[(A -I- g^/2m) 4- c^q^]"^/^. 

To incorporate in the model the presence of randomly 
distributed impurities wc add to Hamiltonian ([!]) a ran- 
dom perturbation which can be parametrized in full gen- 
erality by the four disorder potentials Vi{r): 



H = Ho+ J2 '^'^'W 



(2) 



i—Q.x.y.z 



The potentials Vi{r) are assumed to be random Gaussian 
fields with zero mean Vi{r) = 0. The disorder strengths 
ai parametrize the variance of the distribution according 
to 



%<5(r - r'). 



All these terms have simple physical meanings: Vb(r) is a 
random scalar potential, T4 (r) describes local fluctuations 
of A, Vy{r) is a random gauge potential acting in the y - 
directions, while t4(r) plays the role of the random mass. 
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Fig. 2: The averaged over disorder Green function in the 
SCBA. 

It is worth to note that the scmi-mctal-insulator transi- 
tion in the presence of disorder resembles to some extent 
the critical behavior of a classical system at the Lifshitz 
point [30] . 

The self-consistent Born approximation . To 

describe the effect of disorder on the DOS at low energy, 
we naturally employ the self consistent Born approxima- 
tion (SCBA). The SCBA is expected to be reliable, ex- 
cept possibly in the very vicinity of the transition [24l[29] . 
The retarded and advanced Green functions for the Hamil- 
tonian without disorder ([T]) written in Fourier space are 
given by 



(A + £r 



c^x + cqyay 



(4) 



To calculate the corresponding Green functions averaged 
over disorder within the SCBA one takes into account the 
one-loop diagram contributing to the self-energy. Replac- 
ing the bare Green function in this diagram by the dressed 
one, we obtain a self-consistent equation for the self-energy 
S(£) graphically depicted in figure[2j The equation is sim- 
plified by neglecting the external momenta depcndancc of 
the self-energy S(e). For convenience we introduce the 
matrix function X{e) ~ e ~ ^(e). Making an ansatz 
X{e) = A'o(e)cro — X.\{e)axi where Xi{e) are scalar func- 
tions, we can rewrite the SCBA equation as 

27rca /■ Xo(e) 



^o(e) 



2TTca 



n[x,,qx] 



g n[Xi,qx] 



c^ql 



(5a) 



(5b) 



Here we have denoted := / (P'q/{2nY and r2[Xi,ga;] := 
A 4- Xi{e) + ql/{2m,), and also defined the total disorder 
strength a := ao +C(x +ciy +ciz and the disorder asymme- 
try a := ao + ax — cxy — a^. We first notice that the UV 
diverging integral in eq. (|5b[) is determined by the cutoff 
A. To simplify calculations we impose the cutoff only on 
wavevcctors in the x-direction, so that \qx\ < A, while we 
will neglect it in the y-direction. Integrating over q^ and 
qy in eq. (|5bp we obtain Xi ~ —aA/y/m for large A which 
docs not depend on e. It is convenient to introduce the 
renormalized driving parameter 



A = A- 



A, 



(6) 



which accounts for the shift due to disorder asymmetry of 
the scmi-mctal to insulator transition point from A = 
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Fig. 3: The DOS in units of \/rn/c deep in tlie insulating 
phase for A = 1 (bottom plot) and in the semi-metal phase for 
A = — 1 (upper plot). The disorder strength is a = 0.02. In 
the insulating phase the DOS vanishes for e < £c ~ 0.91. The 
red solid lines are obtained using numerical solution of eqs. (|8} 
and (|12p ; blue dashed lines using the first iteration of the latter 
equations that gives the shifted DOS of the clean system; black 
dots using numerical solution of the simplified eq. ([9]). The 
black dot-dashed line is the power-law approximation (llip . 



towards A = 0. The shift depends on the cutoff and thus 
is non universal similar to what happens at the classical 
Lifshits point [3D]. Solving eq. (|5ap for we find the 
DOS from the expression 



P{e) 



1_ 



-Im Tr / G 



-q) 



Mo«(e). (7) 



The form of the integral over q in eq. ([5a|) strongly depends 
on the sign of A, i.e. on the direction from which we 
approach the topological transition. Hence we first discuss 
the effect of disorder deep in the insulating (A > 0) and 
semi-metal (A < 0) phases, and finally we concentrate on 
the vicinity of the topological transition (A « 0). 

The insulating phase. For A > 0, integrating over q 
in eq. ([5a|) we obtain 



2A 



^0 



(8) 



where K{z) = dt[{l - t^){l - zt"^)]-^/"^ is the com- 
plete elliptic integral of the first kind0 Equation ([5]) 
has two complex conjugate solutions corresponding to the 
retarded and advanced functions. We always find for 
A > Ac = the DOS of an insulator: there ex- 

ists an energy gap ec = A 4- 0(a) such that the solution 
of eq. ([8|) is purely real for < e < Ec, and thus, the 
DOS p{e) vanishes. The DOS computed using eq. ([7]) and 
the imaginary part of the numerical solution of eq. ([8]) is 
shown in figure [3] For \e — ed ^ one can solve eq. ([HI) 
by iterations: this provides an expansion of the solution 
in powers of disorder strength. The first iteration derived 
by taking Xq e + iO+ in the r.h.s of eq. ([5]) gives a lin- 
ear in disorder correction and thus the DOS of the clean 



^ Note that this definition differs from that used in |13) by chang- 
ing K{z) ^ K{y^). 



system found in [T3] up to the shift A ^ A. It is also 
shown in figure [3] for comparison. To study the solution 
of the SCBA equation in the vicinity of the energy gap £c, 
we replace the complete elliptic integral in eq. ([S]) by its 
asymptotics that yields 



Xn 



In 



Xo-^ 



(9) 



where "+/~" correspond to the retarded and advanced 
functions, respectively. Equation ([9]) can be written as 
Xq = e -I- fiXo) so that the energy gap Sc and the corre- 
sponding value Xo{ec) — Xc satisfy: Xc ~ Sc + f{Xc) and 
f'{Xc) = 1. In the limit of weak disorder we obtain 



= A 1 - 



2VA 



1 - In 



64VA 



+ 0{a^lna), (10) 



and the expansion of the DOS around this point for s ^ Sc' 



p{e) 



f'iXc) 



+ O {{e - e.r/') . (11) 



Thus, the DOS of the disordered system in the insulating 
phase vanishes as a power law with the exponent i in 
contrast to jump in the DOS of the clean system (see the 
bottom plot in figure [3|) . 

The semi-metallic phase. In the case A < 0, the inte- 
gration over q in eq. ([5a)) leads to the equation 



Xn^e- 



V2aXo 
V^ + Xo 



K 



2Xn 



A + Xn 



2iK 



A + Xc 



(12) 



The imaginary part of the numerical solution of eq. (|12p 
gives the DOS which is shown in the upper plot of figure [3| 
Iterating eq. (|12P one can obtain the solution for |e— A| ^ 
0. The DOS of the clean system found in [T3] (up to the 
shift ([H])) is recovered by taking the limit Xq — ^ e -|- 10+ 
in the r.h.s of eq. ([12]) and substituting the solution to 
eq. ([7]). This DOS diverges logarithmically at e = |A| 
while the DOS of the disordered system is continuous and 
smooth (sec the upper plot in figure [3]). 

Topological transition. Let us now focus on the region 
of small A. The DOS computed numerically using the 
SCBA eqs. ([8]) and ([12]) is shown in figure [4] as a function 
of A in the vicinity of the transition A = 0. One can see 
that in the presence of disorder there is a region around 
the transition point A = characterized by a finite DOS 
at zero energy, p{£ = 0) > 0. The DOS at e = strictly 
vanishes for A > Ac = tt'^o^ /2. Thus, in the presence 
of disorder the transition (semi-Dirac) point between the 
Dirac semi-metal and the insulator is replaced by a new in- 
termediate regime which we call the disordered semi-Dirac 
regime. The disordered Dirac phase smoothly evolves into 
this intermediate regime, while the SCBA suggests a sharp 
continuous transition at A = Ac between this disordered 
semi-Dirac phase and the insulating phase. Exactly at the 
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Fig. 4: The DOS in units of ^/m/c close to the transition as a 
function of A (upper panel at arbitrary energy, bottom panel 
at zero energy). The disorder strength is a = 0.02. The red 
dashed hue is the DOS at A = given by eq. gS]). The DOS 
at the critical point Ac ~ 0.002 is given by eq. (|13|) . Dotted 
Une in the bottom panel is given by eq. HJ). 



transition wc find an algebraic behavior of the DOS 



1/6 



(7ra) 



l/3pl/3 



(13) 



which is displayed in figure ID The extrapolation towards 
the disordered semi-metallic regime is obtained for large 
negative A < 0: in this regime we recover a DOS at zero 
energy decaying as 



p(0) = 



lAI 



"^g-VjA|/(av^)^ 



(14) 



Hence, instead of a transition the system undergoes a 
crossover to disordered Dirac fermions. Indeed, this be- 
havior is fully consistent with the SCBA picture for the 
disordered Dirac fermions. For graphene it predicts a 
finite DOS p[e) = 2Tq/ {■K'^vlao) for e < Tq where 



To = A^ie 



is the energy scale generated by disor- 



der. Here A^i is the bandwidth, - the Fermi velocity 
and ao - the dimensionless strength of disorder [24] . This 
would imply a finite DOS at zero energy for all considered 
types of disorder. However, this contradicts the results ob- 
tained for random gauge and random mass disorder using 
more reliable methods. For instance, for random gauge 
disorder one expects p{e) ~ e^/^~^ with z = 1 
weak disorder case {ao < 2) [20] and z = [Sao] 
in strong disorder case {an > 2) [22]. Similarly, with ran- 
dom mass disorder, the DOS remains linear in e up to 
logarithmic corrections [26]. Only for scalar potential dis- 
order the DOS saturates at a finite value in the vicinity of 
the Dirac point [24] . It was argued in [21] that the failure 
of SCBA in the vicinity of the Dirac point for e < Fo is 



m 

^1/2 _ 1 



due to importance of the diagrams with crossed disorder 
lines, neglected within the SCBA which takes into account 
only the non-crossed ones. Another possible reason for 
the failure of the SCBA is the divergence of the fermion 
wavelengths at the Dirac point rendering the SCBA un- 
controllable, i.e. without a small parameter. Nevertheless 
one can still rely on the SCBA for energies e ^ Fq. The 
question of reliability of the SCBA in the vicinity of the 
topological transition is even more subtle. 

Disordered Semi-Dirac. It is interesting to note that 
the SCBA equations can be analytically solved exactly 
at the point A = 0, inside the "disordered semi-Dirac" 
regime. Indeed using that K{2 + 0~^/z) — C(l — 
sign[Im z ]i)/V2 where C = T{l/Af/{4y^) one can 
rewrite eq. (jH]) as ATq = e + (1 T i)Ca^/-XQ , where "-/-H" 
corresponds to the retarded/advanced Green's function. 
We obtain the DOS at the point A = 



Pis) 




(15) 



such that p(0) = 2y/mC'^a/ {ctt'^). Note that the expan- 
sion of cq. (|T5]) in powers of the disorder strength a is in 
fact expansion in aj^fe. Its truncation to a finite order 
in a appears to be pathological for e <C a^. We interpret 
this as a sign of the probable breakdown of the SCBA 
for a small enough energy. In this regime, a more refined 
technique is required [25] . 

Before we consider a weak disorder renormalization 
group (RG) approach to the topological transition let us 
first discuss the consequences of our SCBA results for 
the behavior of specific heat and magnetic susceptibil- 
ity [13l[3T]. Assuming that the chemical potential is tem- 
perature independent and fixed at zero energy the specific 
heat can be written as eq. (18) in [T3]. In the clean system 
the low T specific heat interpolates from a behavior 
far in the semi-metal phase to an activated behavior in 
the insulating phase. Exactly at the transition it exhibits 
a T^l"^ behavior. In the presence of disorder the low T 
specific heat at the critical point should behave as T^l'^ . 
The magnetic susceptibility x is related to the DOS by 
eq. (10.11) in [32]. In the clean system, x is activated in 
the insulator, x ^ 7" in the semi-metal phase and x ^ T^/^ 
at the transition. With disorder, the scaling behavior of 
the susceptibility at the transition is modified to x ~ T^l'^ . 

Weak disorder renormalization group. We now 

study whether the critical point found using the SCBA is 
accessible within a weak disorder RG approach. The start- 
ing point is the generating functional for the Green func- 
tions Z = J Dil;Dil;e~^ , where S is the effective action 
at a given energy (Matsubara frequency) e. The effective 
action can be written in terms of two fcrmionic (anticom- 
muting) fields ip{r) and ip{r) which are two completely 
independent Grassman variables. For generality we con- 
sider a class of d dimensional models in which the bare 
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dispersion is quadratic in I direction and linear in d — I 
remaining directions. For instance, the case {d = 2,1 = 1) 
describes the semi-Dirac fermions, {d = 3,1 = 0) - Weyl 
fermions, and {d = 3,1 = 1) - semi- Weyl fermions [55] . 

Scaling dimensions. The standard way to general- 
ize the model to d dimensions is to introduce the gen- 
erators of the d-dimensional Clifford algebra. In even 
dimensions these generators can be represented by the 
matrices ji (i = l,...,d) of dimension 2"^/^. Generally 



one also defines 7^+1 



7s 



71... 7c;. For the 



case of odd dimensions one constructs the matrices 7^ 
(i = l,...,d — 1) similarly to the case of even dimension 



and adds 7^ = 75 



-(d-l)/2. 



7i...7d_i. These matrices 



satisfy the anticommutation relations: 

111] + = 2%7o, «, .7 1, [M + 1]] , 



(16) 



where 70 = I is the identity matrix and [[d + 1]] stands 
for d + 1 if is even, and d if it is odd. The matri- 
ces 7i and 75 replace the Pauli matrices cr^ (i = x, y) 
and (Tz in d dimensions. Wc define the two matrix vec- 
tors 7|| := (71,..., 7;) and 7j_ (7;+!, 7^) and replace 
in Hamiltonian ^ dx ^ 9|[ which acts in /-dimensional 
subspacc and dy — >■ d±^ which acts in [d — /)-dimensional 
subspace. Introducing n copies of the original system and 
using the replica trick to average over disorder we can 
write the effective replicated action as 



s = j d'v\y^i,a 



TTC 

fm 



.a=l 
n [[d+1]] 



711 A - 



2m 



ic7j_(?j_ — ie 



E o.AMrhiMr)][MrhjMr)]}, (17) 



a, 6=1 j=0 



where d r d'-rn d"^ 'r±. The properties of the initial 



disordered model are obtained in the limit n ^ 0. Dimen- 
sional analysis of the action p7)) yields the bare dimen- 
sions: [rii] = i-i^^''^, [r±] = i-L^^, 



and [aj 
and 



jj^ , where /i is an arbitrary momentum scale 



S = 2 + l/2-d. 



(18) 



Thus the upper critical dimension of the model (jl7p is 
dc = 2 + 1/2 and disorder is (i) relevant for S > 0, e.g. 
for 2D semi-Dirac fermions with d — 1/2; (ii) marginally 
relevant for 2D Dirac fermions with 6 — 0; (iii) irrelevant 
for 3D Weyl {5 = —1) and semi- Weyl fermions {S ~ —1/2). 

One-loop RG with restricted disorder. We focus here 
on the case of 2D semi-Dirac fermions. To renormalize 
the model ([TT]) one has to eliminate UV divergences in 
the loop expansion by introducing counter terms. The 
one-loop counterterms to the propagator and disorder are 
given by the diagrams shown in figure [5l To one- loop or- 
der the m and c do not get corrected while A gets an 
additive shift similar to that in eq. ©. We use dimen- 
sional regularization in which this shift vanishes |30| . The 
diagrams (c) and (d) in figure [5] can generate extra terms 



(a) 



(b) 



(c) 



(d) 



Fig. 5: The one-loop diagrams contributing to the propagator 
(a) and disorder renormalization (b)-(d). 



with new algebraic structures like 7i7j7fe ® Jijjjk and so 
on. In d = dc — 6 dimensions they have to be treated 
as independent new interactions reflecting the fact that 
within dimensional regularization this theory is not mul- 
tiplicatively renormalizable |34j . For the Dirac fermions 
this problem shows up to orders higher than one-loop and 
can be solved, for example, using a cutoff in strictly 2D so 
that the extra interactions do not arise. In contrast, for 
the semi-Dirac fermions the extra interactions infinitely 
proliferate in the renormalization already to the one-loop 
order. This cannot be cured by using a different regular- 
ization scheme since the upper critical dimension | is non 
integer. To avoid this obstacle we restrict ourselves to the 
presence of a single type of disorder. It easy to see that 
the diagrams (c) and (d) in figure [5] with any type of dis- 
order always contribute to ai = a^. Thus we can derive a 
closed one-loop fiow equation by considering a model with 
only the ax type of disorder. The flows of the renormal- 
izcd dimensionless energy and disorder strength read for 
this model 



—f-id^ Ine = 1 
—^df^ax = S a 



(19a) 
(19b) 



where we have included the value of the one-loop inte- 
gral in redefinition of ax- Equation ()19bl) possesses a 
fixed point a* = 5. Integrating eq. ()19ap up to the 
scale corresponding to e ~ 1, wc find e ~ 11^ with 
z = l + a*x = l + 5+0{d^). The DOS behaves as p - {-^jip). 



which leads to its scaling as p ■ 
terms of energy as 



IJ-' 



3/2- 



or equivalently in 



P(e) 



.(3/2-z)/2 



,(l/2-5)/(l+5)+0(5^) 



(20) 



Hence we find that for 5 = ^ the DOS exponent vanishes to 
one loop order suggesting a constant DOS at zero energy 
when only the ax disorder is present. This result improves 
upon the critical behavior extracted from the SCBA anal- 
ysis which was shown to be invalid at low energy. Whether 
the exponent vanishing is accidental or holds beyond the 
one-loop order remains an open question. 

Full one-loop RG analysis. To analyze the scaling be- 
havior of the full disordered 2D semi-Dirac fermions, we 
find it more technically tractable to consider an alternative 
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generalized model given by the effective action 



S = j d^'vH^i^a 



n 

^ ^ ^ aj[?/;a(r)crj?/'a(r)][V^6(r)crj^/.6(r)]|. 



a, 6—1 j— 0,a:,y,2: 

where (fir := d''~-'^rj| drj^. The advantage of this model is 
that it involves no infinite Clifford algebra. Dimensional 
analysis of this model suggests that its upper critical di- 
mension is dc = 3. We can construct a perturbative in 
disorder RG around this dimension using expansion in a 
small parameter 5 = (3 — (i)/2 similar to that done for 
the model ([TT]). The corrections to energy and disorder 
strengths are given by the same diagrams shown in fig- 
ure The one- loop RG equations read 

—^id^ h\e ^ I + + + cty + az, 

-^idf^UQ (5q;o + 4q;o(q!o + Zofj; + ay + a^) + ^ia^az, 



-^d^ay ^ Say 
—^d^az ~ 6az — 4:az{ao — 2ax — ay 



Aaoax- 



Unfortunately, the solutions of these RG equations always 
flow towards a strong disorder regime which describes An- 
derson localization, but is not within the reach of a weak 
disorder RG analysis. In this generalized model, we do no 
flnd a critical perturbative fixed point, as opposed to the 
model with only ax disorder. 

Conclusions. — We have studied the effect of disor- 
der on the topological merging transition from a Dirac 
semi-metal to a band insulator in 2D systems. Using the 
SCBA we have found that disorder smears out the tran- 
sition and an intermediate disordered semi-Dirac regime 
emerges. This phase is separated from the insulating phase 
by a sharp continuous transition while there is a crossover 
between this phase and disordered Dirac fermions regime. 
Exactly at the transition the SCBA suggests a power-law 
behavior of DOS with a nontrivial exponent i. We have 
also shown that the critical point can be perturbatively 
described using a weak disorder RG when only random 
fluctuations of A are present but is unfortunately inacces- 
sible for the most general model. 
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